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Abstract 

Despite the theory of neutrino oscillations being rather old, some of its basic issues are 
still being debated in the literature. We discuss, in the framework of the wave packet 
approach, a number of such issues, including the relevance of the "same energy" and 
"same momentum" assumptions, the role of quantum-mechanical uncertainty relations 
in neutrino oscillations, the dependence of the production/detection and propagation 
coherence conditions that ensure the observability of neutrino oscillations on neutrino 
energy and momentum uncertainties, the question of (in)dependence of the oscillation 
probabilities on the neutrino production and detection processes, the applicability lim- 
its of the stationary source approximation, and Lorentz invariance of the oscillation 
probability. We also develop a novel approach to calculation of the oscillation prob- 
ability in the wave packet picture, based on the summation/integration conventions 
different from the standard one, which gives a new insight into the oscillation phe- 
nomenology. We discuss a number of apparently paradoxical features of the theory of 
neutrino oscillations. 
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1 Introduction 



More than 50 years have aheady passed since the idea of neutrino oscillations was put 
forward [1, 2], and over 10 years have passed since the experimental discovery of this phe- 
nomenon [3]. However, surprisingly enough, a number of basic issues of the theory of 
neutrino oscillations are still being debated. Moreover, some features of the theory appear 
rather paradoxical. The issues that are still under discussion include 

(1) Why do the often used same energy and same momentum assumptions for neutrino 
mass eigenstates composing a given flavour state, which are known to be both wrong, 
lead to the correct result for the oscillation probability? 

(2) What is the role of quantum-mechanical uncertainty relations in neutrino oscillations? 

(3) What determines the size of the neutrino wave packets? 

(4) How do the neutrino production/detection and propagation coherence conditions that 
ensure the observability of neutrino oscillations depend on neutrino energy and mo- 
mentum uncertainties? 

(5) Are wave packets actually necessary for a consistent description of neutrino oscilla- 
tions? 

(6) When can the oscillations be described by a universal (i.e., production and detection 
process independent) probability? 

(7) When is the stationary source approximation valid? 

(8) Would recoillessly emitted and absorbed neutrinos (produced and detected in Mossbauer- 
type experiments) oscillate? 

(9) Are oscillations of charged leptons possible? 

In the present paper we consider the first seven issues listed above, trying to look at them 
from different perspectives. We hope that our discussion will help clarify these points. For 
the last two issues, we refer the reader to the recent discussions in [4, 5, 6] (for oscillations of 
Mossbauer neutrinos) and [7] (for oscillations of charged leptons). We also discuss Lorentz 
invariance of the oscillation probability in the wave packet picture of neutrino oscillations. 

The structure of the paper is as follows. In Sec. 2 we review the standard derivations of 
the neutrino oscillation probability in the plane wave and stationary state approaches and 
point out their inconsistencies. In Sec. 3 we introduce a general wave packet formalism for 
neutrino oscillations which does not rely on the specific form of the wave packets and discuss, 
on its basis, why the "same energy" and "same momentum" approaches give the correct 
oscillation probability. The reader not interested in the latter question can skip Sec. 3.3. 
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In Sec. 4 we develop a novel approach to calculation of the oscillation probability in the 
wave packet picture, based on the alternative summation/integration conventions. Sec. 5 
is devoted to the discussion of the role of the quantum-mechanical uncertainty relations in 
neutrino oscillations. Here we consider, in particular, the question of what determines the 
size of the neutrino wave packets, as well as the conditions for coherent neutrino production 
and detection and the issues of coherence loss on the way between the neutrino source and 
detector and its possible restoration at detection. In Sec. 6 we consider the stationary source 
approximation for neutrino oscillations; we derive this approximation in the wave packet 
picture and analyze its domain of applicability. In Sec. 7 we discuss the question of when the 
oscillation probability is independent of the neutrino production and detection processes. 
Sec. 8 is devoted to the discussion of the Lorentz invariance of the oscillation probability. 
We add a few comments on the wave packet approach in Sec. 9 and summarize our results 
in Sec. 10. Appendices A and B contain some auxiliary material. Sections 4, 6, 7 and 8 are 
largely independent of each other and can be read separately. 

We plan to update this article in the future if and when new unresolved questions or 
paradoxes of the theory of neutrino oscillations come to light. 

2 Same energy or same momentum? 

In most derivations of the so-called standard formula for the probability of neutrino oscilla- 
tions in vacuum (see eq. (6) below), usually the assumptions that the neutrino mass eigen- 
states composing a given flavour eigenstate either have the same momentum [8, 9, 10, 11] 
or the same energy [12, 13, 14, 15] are made. The derivation typically proceeds as follows. 

First, recall that in the basis in which the mass matrix of charged leptons has been 
diagonalized the fields describing the massive neutrinos Vi and flavour-eigenstate neutrinos 
Pa and the corresponding states and |z/a) are related by 



where U is the leptonic mixing matrix. If one now assumes that all the mass eigenstates 
composing the initially produced flavour state |i^(0)) = \va) have the same momentum, 
then, after time t has elapsed, the ith mass eigenstate will simply pick up the phase factor 
exp(— ■iE'jt), and the evolved state \v{t)) will be given by 



Projecting this state onto the flavour state \ub) and taking the squared modulus of the 
resulting transition amplitude, one gets the probability of the neutrino flavour transition 
after the time interval t: 





(1) 




(2) 




(3) 
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Next, taking into account that the energy Ei of a relativistic neutrino of mass and 
momentum p is 



2 

mf 



Ei = ^Jp^ + mf ^ (4) 

and that for relativistic pointhke particles the distance L they propagate during the time 
interval t satisfies 

L ^ t, (5) 

one finally finds 

I Am?. 2 

P{va^Vb^L) = \Y,Uue-'^''U:, , (6) 



where Arnfj = mf — and the index j corresponds to any of the mass eigenstates. This is 
the standard formula describing neutrino oscillations in vacuum.^ Note that in this approach 
neutrino states actually evolve only in time (see eq. (3)); the usual coordinate dependence 
of the oscillation probability (6) is only obtained by invoking the additional "time-to-space 
conversion" assumption (5). Without this conversion, one would have come to a paradoxical 
conclusion that neutrino oscillations could be observed by just putting the neutrino detector 
immediately next to the source and waiting long enough. 

Likewise, one could assume that all the mass-eigenstate neutrinos composing the initially 
produced flavour state |z/a) have the same energy. Using the fact that the spatial propagation 
of the ith mass eigenstate is described by the phase factor e'^*^ and that for a relativistic 
neutrino of mass rrii and energy E 



p. = ^E^ -mfc^E- mf/2E , (7) 

one again comes to the same standard formula (6) for the oscillation probability. Note that 
in this case the neutrino flavour evolution occurs in space and it is not necessary to invoke 
the "time-to-space conversion" relation (5) to obtain the standard oscillation formula. 

The above two alternative derivations of the oscillation probability are very simple and 
transparent, and they allow one to arrive very quickly at the desired result. The trouble 
with them is that they are both wrong. 

In general, there is no reason whatsoever to assume that different mass eigenstates 
composing a flavour neutrino state emitted or absorbed in a weak-interaction process have 
either the same energy or the same momentum ^. Indeed, the energies and momenta of 
particles emitted in any process are dictated by the kinematics of the process and by the 
experimental conditions. Direct analysis of, e.g., 2-body decays with simple kinematics, such 
as TT^ —> /^ + z//(z/;), where I = e, fi, allows one to find the 4-momenta of the emitted particles 



^Let us stress that in this paper by the standard oscillation probability we will always understand the 
non-averaged (i.e., p- and L-dependent) probability (6). 

^The only exception we are aware of are neutrinos produced or detected in hypothetical Mossbauer-type 
experiments, since for them the "same energy" assumption is indeed justified very well. 
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and shows that neither energies nor momenta of the different neutrino mass eigenstates 
composing the flavour state ui are the same [16, 17]. One might question this argument 
on the basis that it relies on the energy-momentum conservation and the assumption that 
the energies and momenta of the emitted mass eigenstates have well defined (sharp) values, 
whereas in reality these quantities have intrinsic quantum-mechanical uncertainties (see the 
discussion in Sec. 5.1). However, the inexactness of the neutrino energies and momenta does 
not invalidate our argument that the "same energy" and "same momentum" assumptions 
are unjustified, and in fact only strengthens it. It should be also noted that the "same 
energy" assumption actually contradicts Lorentz invariance: even if it were satisfied in 
some reference frame (which is possible for two neutrino mass eigenstates, but not in the 
3-species case), it would be violated in different Lorentz frames [17, 18]. The same applies 
to the "same momentum" assumption. 

One may naturally wonder why two completely different and wrong assumptions ( "same 
energy" and "same momentum" ) lead to exactly the same and correct result - the standard 
oscillation formula. To understand that, it is necessary to consider the wave packet picture 
of neutrino oscillations. This will also allow us to analyze many other issues of the theory 
of neutrino oscillations. 

3 Shape-independent wave packet approach to neu- 
trino oscillations 

3.1 Neutrino wave packets and the oscillation probability 

In the discussion in the previous section we were actually considering neutrinos as plane 
waves or stationary states; strictly speaking, this description was inconsistent because such 
states are in fact non-propagating. Indeed, the probability of finding a particle described by 
a plane wave does not depend on the coordinate, while for stationary states this probability 
does not depend on time. In quantum theory propagating particles must be described 
by moving wave packets (see, e.g., [19]). We will therefore consider now the wave packet 
approach to neutrino oscillations. Many features of this approach have been discussed 
previously [20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32]. However, unlike in the majority 
of those papers, where Gaussian wave packets were considered, in our analysis we will not 
assume any specific shape of the wave packets. 

Let a flavour eigenstate Ua be produced during a time interval Ats centered at t = 
in a source centered at x = 0. The wave packet describing the evolved neutrino state at a 
point with the coordinates {t, x) is then 



We consider the evolution of the neutrino state (8) after the formation of the wave packet 




(8) 
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has been completed, that is at t > Ats- This state carries an information on the neutrino 
production process and corresponds to the mass-eigenstate components of the flavour neu- 
trino state being on the mass shell (see Sec. 5.2 for the discussion of this point). In what 
follows, we shall call it the production state. In eq. (8), \E'i(x, t) is the wave packet describing 
a free propagating neutrino of mass rrii: 



^i(x, t) 



(27r)3/2 



/f (P - P.) e*''""''''^^^* 



(9) 



where ff{p — Pi) is the momentum distribution function with pj being the mean momentum, 
and Ei{p) = \/p'^ + mf. The superscript S at /f (p — Pi) indicates that the wave packet 
corresponds to the neutrino produced in the source. We will assume the function /f (p — Pi) 
to be sharply peaked at or very close to zero of its argument (p = p^), with the width of 
the peak aps <^ Pi- ^ No further properties of /f (p — Pi) need to be specified. Expanding 
E.i{p) around the mean momentum. 



Ei{p) = EM 



dE,{p) 



dpi 



{p - piY + 



Id'Eiip) 



2 dpWp'' 



{p - p,y {p - p,f + ... 



(10) 



one can rewrite eq. (9) as 



v^,(x, t)^e^P^''-^^'(^')*^7f(x-v,,t), 



'11^ 



where 



is the shape factor and 



t) 



S'p 



(27r)3/2 



/f (p) e^P(^-^«'*) 



dE, 



' gi 



dp 



E,, 



(12) 



(13) 



is the group velocity of the wave packet. Here we have retained only the first and the second 
terms in the expansion (10), since the higher order terms are of the second and higher order 
in the small neutrino mass and so can be safely neglected in all the situations of interest, 
with a possible exception of supernova neutrinos. This approximation actually preserves 
the shape of the wave packets (and, in particular, neglects their spread). Indeed, the shape 
factor (12) depends on time and coordinate only through the combination (x — Vgit); this 
means that the wave packet propagates with the velocity v^j without changing its shape. 

If the momentum dispersion corresponding to the momentum distribution function 
/f (p — Pi) is aps, then, according to the Heisenberg uncertainty relation, the length of 



the wave packet in the coordinate space axs satisfies o"; 



xS 



the shape factor function 



•^For symmetric wave packets (i.e. when ff{p — Pi) is an even function of its argument), the position 
of the center of the peak coincides with the mean momentum p^. For asymmetric wave packets, it may be 
displaced from p^. 
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(^f (x — Vgjt) decreases rapidly when |x — Vgjt| exceeds (Txs- Eq. (11) actually justifies and 
corrects the plane-wave approach: the wave function of a propagating mass-eigenstate neu- 
trino is described by the plane wave corresponding to the mean momentum pj, multiplied 
by the shape function f^f (x — v^jt) which makes sure that the wave is strongly suppressed 
outside a finite space-time region of width ax around the point x = Wgit. 

In the approximation where the wave packet spread is neglected, the evolved neutrino 
state is given by eq. (8) with \&j(x, t) from eq. (11). Note that the wave packets corresponding 
to different neutrino mass eigenstates Vi are in general described by different momentum 
distribution functions /j^(p — Pi) and therefore by different shape factors fiff (x — v^jt). 

Let us now turn to the detected flavour-eigenstate neutrino z/f,. We describe its state by 
a wave packet peaked at the coordinate L of the detecting particle: 

|z.,(x-L)) = ^f/*v[/f(x-L)|z.,). (14) 

i 

This state has no time dependence because the detection process is essentially time inde- 
pendent on the time scale of the inverse energy resolution of the detector (see Sec. 6). The 
wave function (x — L) can be written as 

(X - L) = / //'(P - P^) e^''^''"'^ > (15) 

where //^(p — p^) is the momentum distribution function of the wave packet characterizing 
the detection state, with p- being the mean momentum. Note that in general the average 
momenta of the production and detection states corresponding to the same neutrino mass 
eigenstate, pj and p-, are different. This is because the detection process may be sensitive to 
the interval of neutrino energies which does not exactly coincide with the energy spectrum 
of the emitted neutrino state (we shall discuss this point in more detail below). Just like 
for the production state, we will assume the function //^(p — p-) to be sharply peaked at 
or very close to zero of its argument (p = p^), with the width of the peak (jpD ^ p'i- By 
shifting the integration variable in eq. (15) one can rewrite ^E'f (x — L) as 

*f(x-L)=e^P^(''-^)<7f(x-L), (16) 

with 

jf(--L) = /(^ (17) 

The function fiff (x — L) is the shape factor of the wave packet corresponding to the detection 
of the 2th mass eigenstate. It decreases quickly for |x — L| > axu, where Gxd is the spatial 
length of the detection wave packet which is related to the momentum dispersion by 
the Heisenberg relation Gxd ^ 
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The transition amplitude Aabi^it) describing neutrino oscillations is obtained by pro- 
jecting the evolved state (8) onto (14): 



Aab{'L,t) = j d'x{Ub{^-l.)\Ua{^,t)) = Y,U:rUu ^ rf'x * (x - L) (x, t) . (18) 

Substituting here eqs. (11) and (16) yields 

A6(L, t) = Y. ^-^^^ ^^(L - v,.^) e-^^»(^'»)*+*P'^ , (19) 

i 

where 

^.(L - wg^t) = j d'xgfi^ - v,,t)(7f *(x - L) e^(p.-P^)(-L) . (20) 

That this integral indeed depends on L and Vgit only through the combination (L — v^jt) 
can be easily shown by shifting the integration variable in (20). Gj(L — v^jt) is an effective 
shape factor whose width cr^ depends on the widths of both the production and detection 
wave packets axs and a^D and is of the order of the largest between them. Indeed, since 
the moduli of the shape factor functions ' quickly decrease when the arguments of these 
functions exceed the corresponding wave packet widths a^s or a^D, from eq. (20) it follows 
that Gi(L — Vgit) decreases when |L — Vgjt| becomes large compared to max{axs, ^xd}- 
Actually, since characterizes the overlap of the wave packets describing the production 
and detection states, it exceeds both a^s and ct^d. In particular, for Gaussian and Lorentzian 
(in the coordinate space) wave packets one has = \/ + a^j^ and ax = (Jxs + o'xD, 
respectively. 

The probability Pafe(L, t) = P{va i^b] L, t) of finding a flavour eigenstate neutrino Uf, 
at the detector site at the time t is given by the squared modulus of the amplitude ^^(L, t) 
defined in eq. (18). Since in most experiments the neutrino emission and arrival times are 
not measured, the standard procedure in the wave packet approach to neutrino oscillations 
is then to integrate Pab(L, t) over time. In doing so one has to introduce a normalization 
factor which is usually not calculated, ^ and in fact is determined by imposing "by hand" the 
requirement that the probabilities Pab{L) satisfy the unitarity condition. This is an ad hoc 
procedure which is not entirely consistent; the proper treatment would require to consider 
the temporal response function of the detector and would automatically lead to the correct 
normalization of the oscillation probabilities. It is shown in Sec. 6 that taking into account 
the temporal response function of the detector confirms the standard result, and therefore 
we follow the same procedure here. The proper normalization of the oscillation probability 
is achieved by imposing the normalization condition 

/oo 
rft|G',(L-v,,t)|2 = l. (21) 
-oo 

^For an exception, see [2G]. 
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For simplicity, from now on we neglect the transverse components of the neutrino mo- 
mentum, i.e. the components orthogonal to the line connecting the centers of the neutrino 
source and detector; this is a very good approximation for neutrinos propagating macro- 
scopic distances. The probability of finding a z/f, at the detector site provided that a Ua was 
emitted by the source at the distance L from the detector is then 

/oo 
dt\AabiL,t)\^ = Y,U:^brUakU;jik{L) , (22) 

i,k 

where 

/oo 
dtG,{L - vg,t) Gl{L - Vgkt) e-*^'^-^^'*) . (23) 
'OO 

Here Gi{L — Vgit) is the effective shape factor corresponding to the zth neutrino mass eigen- 
state defined in (the 1-dimensional version of) eq. (20). The quantity A(f)ik{L, t) is the phase 
differences between the ith and kth mass eigenstates: 

A(j),k = {Ei - Ek)t - {pi - pk)L = AEik t - Apik L , (24) 

where 

E, = ^p^ + ml (25) 

Note that this phase difference is Lorentz invariant. 

To calculate the observable quantities - the numbers of the neutrino detection events 
- one has to integrate the oscillation probability (folded with the energy spectrum of the 
source, detection cross section and detector efficiency and energy resolution functions) over 
the neutrino spectrum and over the macroscopic volumes of the neutrino source and detector. 



3.2 The standard oscillation formula and beyond 

According to eq. (22), the oscillation probability can be expressed through the elements 
of the leptonic mixing matrix and the integral Iik{L) defined in eq. (23). Let us now 
consider the properties of this integral. First, we derive a useful representation for the 
quantities Gi{L — Vgit). Expressing the shape factors g^'^{x) of the wave packets through the 
corresponding momentum distribution functions according to (the 1-dimensional versions 
of) eqs. (12) and (17) and substituting the result into (20), we find 

/OO 
dpf!{p)fF*{P + ^^)e'''^''-''^''\ (26) 
'OO 

where 

Si=Pi- Pi (27) 

is the difference of the mean momenta of the production and detection states. If 6i exceeds 
significantly the sum of the widths of the momentum distributions of the production and 
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detection states aps+CTpD, the integral in (26) is strongly suppressed due to the lack of overlap 
of the functions ff{p) and fP*{p + Si) in the integrand. This is merely a manifestation 
of the approximate momentum conservation in the process: the mean momentum of the 
detected neutrino state should be approximately equal to that of the emitted neutrino, with 
possible deviation not exceeding the overall momentum uncertainty of the production and 
detection processes. Note that in the plane wave limit the momentum distribution functions 
fiip — Pi) and ffip — p'i) are reduced to the corresponding (5-functions, and the momentum 
conservation becomes exact. In what follows we will always be assuming that 6i satisfies 

\Si\ < Crp5 + apD- 

Consider now the integral Iik{L). Substituting eqs. (26) and (24) into eq. (23) we find 

/C30 /"OO /"OO 

dt dpJ dp2ff{pi)f^*{pi+6,)ftiP2)h''{P2+Si) 
-oo J —oo J —oo 

^^ipi{L-Vgit)-ip2{L-Vgkt)+i{Apit:L-AEikt) ^28) 

Performing first the integration over time and making use of the standard integral repre- 
sentation of Dirac's (5-function, we obtain 

UL) = e^(Ap..-Ai.../.,)L _ / dpf!{p)f^*{p + 5.)/f (rp + AE,,/v,u) 

'^gk J -oo 

xf^irp + AE,k/vgk + Si) e^P(^-^)^ , (29) 

where 

Vg = {vg^ + Vgk)/2, r^^-l, (30) 

Vgk 

and we have neglected terms of order (Am^^)^ in the phase factor in front of the integral. 
It can be readily shown that this factor actually contains the standard oscillation phase 
(Am^^,/2p)L (see eq. (34) below), and thus we finally obtain 

UL) = e-^^^ — / dpf!{p)f^*{p + S,)ft{rp + AE,,/vg,) 

^gk J —oo 

xf^irp + AE,k/vgk + S,) 6*^(1-'-)^ . (31) 

We shall now use this representation to discuss the properties of Iik{L). First, we notice 
that in the limit when the group velocities of the wave packets corresponding to different 
mass eigenstates are exactly equal to each other, r = 1, the integral on the right hand side 
of eq. (31) does not depend on the distance L. Since the dominant contribution to this 
integral comes from the region \p\ < ap = 'min{aps, CTpo} of the integration interval, the 
integral is practically independent of L even for r ^ 1 provided that |1 — r\Lap <C 1, or 

i^</coh = ax-J^, (32) 

AVg 
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where Aw^, = \vgi — Vgk\ and ax = l/cp. This is merely the condition of the absence of 
the wave packet separation: the distance traveled by neutrinos should be smaller than the 
distance over which the wave packets corresponding to different mass eigenstates separate 
due to the difference of their group velocities and cease to overlap. If the condition opposite 
to that in eq. (32) is satisfied, the integral hk^L) is strongly suppressed because of the fast 
oscillations of the factor e*^^^"^-*'^ in the integrand. Thus, Iik{L) indicates how well the wave 
packets corresponding to the ith and kth neutrino mass eigenstates overlap with each other 
upon propagating the distance L from the source. 

In addition, Iik{L) is quenched if the split of the arguments of f^'^ and f^'^ in the 
integrand exceeds cxp; therefore, a necessary condition for unsuppressed Iik{L) is 

lAEikWx/vg^l. (33) 

As we shall show in Sec. 5.3, this condition is related to the coherence properties in the neu- 
trino production and detection processes; therefore we shall call it the production/detection 
coherence condition or, for brevity, the interaction coherence condition. 

From the definition ax = CTp^ = [min{aps, o-p£)}]~^ it follows that this quantity is of 
the order of max{ascs, o'xd}, i-e. of the same order of magnitude as the effective length a^ 
of the wave packets Gi{L — Vgit) introduced in the previous subsection. This, in particular, 
means that one can substitute a^ for ax in the coherence conditions (32) and (33). We shall 
be using this fact in what follows. 

If the propagation and interaction coherence conditions (32) and (33) are satisfied, one 
can set r = 1 as well as neglect the quantity AEik in the integral in eq. (31). As shown 
in Appendix A, the limit r = 1 (i.e. Vgi = Vgk) also imphes f^'^{p) = fk'^{p)- From the 
normalization condition (21) and eq. (26) it then follows directly that the integral in eq. (31) 
(with the factor 2T\jvg included) is simply equal to unity, so that Iik{L) = e"**^^™'*/^^'*'^. 
Substituting this into eq. (22) yields the standard oscillation probability. 

3.3 Oscillation phase: answers to the questions 

Since neutrino oscillations occur due to the increasing phase difference between different 
neutrino mass eigenstates, one can learn a great deal about the oscillation phenomenon by 
studying the oscillation phase. We therefore concentrate on this phase now. This will allow 
us to answer some of the question raised in the Introduction and also to look at the results 
of the previous subsection from a slightly different viewpoint. 

Consider the phase difference A0jfc (24) that enters into the expression (23) for Iik{L). 
To simplify the notation, we suppress the indices i and k where it cannot cause a confusion, 
so that AE = AEik, Am^ = Am^^, etc. We shall consider the case \AE\ <^ E which 
corresponds to relativistic or quasi-degenerate neutrinos. In this case one can expand the 
difference of the energies of two mass eigenstates in the differences of their momenta and 
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masses. Retaining only the leading terms in this expansion, one gets 

dE OE 1 
= 'd^^P + d^^"^' = '''^P^2E ' ^^^^ 

where Vg is the average group velocity of the two mass eigenstates and E is the average 
energy. Substituting this into eq. (24) yields [33, 27, 28] 

A</.= ^t - (L-V)AP- (35) 

Note that our use of the mean group velocity and mean energy of the two mass eigenstates in 
eq. (35) is fully legitimate. Indeed, going beyond this approximation would mean retaining 
terms of the second and higher order in Am^ in the expression for A0. These terms are 
small compared to the leading (9(Am^) terms; moreover, though their contribution to A0 
can become of order one at extremely long distances, the leading contribution to A0 is then 
much greater than one, which means that neutrino oscillations are in the averaging regime 
and the precise value of the oscillation phase is irrelevant. 

Let us now consider the expression (35) for the phase difference A0. If one adopts 
the same momentum assumption for the mean momenta of the wave packets representing 
the different mass eigenstates, Ap = 0, the second term on the right hand side disappears, 
which leads to the standard oscillation phase in the "evolution in time" picture. If, in 
addition, one assumes the "time-to-space conversion" relation (5), the standard formula for 
the L-dependent oscillation phase is obtained. 

Alternatively, instead of expanding the energy difference of two mass eigenstates in the 
differences of their momenta and masses, one can expand the momentum difference of these 
states in the differences of their energies and masses: 

Ap = %AE + -^Am" = -AE-^Am\ (36) 
oE am^ Vg 2p 

where p is the average momentum. Substituting this into eq. (24) yields [34] 

Am^ 1 

A0 = — — L (L - Vgt)AE . (37) 

Zp Vg 

Note that this relation could also be obtained directly from eq. (35) by making use of 
eq. (34). If one now adopts the same energy assumption for the mean energies of the wave 
packets, AE = 0, the second term on the right hand side vanishes, and one arrives at the 
standard oscillation phase. 

However, as we shall show now, eqs. (35) and (37) actually lead to the standard os- 
cillation phase even without the same energy or same momentum assumptions. For this 
purpose, let us generically write eqs. (35) and (37) in the form 

A0 = A0,j + A0' , (38) 
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where A(t)st is the standard oscillation phase either in "evolution in time" or in "evolution 
in space" approach, and A0' is the additional term (the second term in eq. (35) or (37)). 
The first thing to notice is that A0' vanishes not only when Ap = (in eq. (35)) or AE = 
(in eq. 37)), but also at the center of the wave packet, where L = Vgt. Away from the 
center, the quantity L — Vgt does not vanish, but it never exceeds substantially the length of 
the wave packet ax-, since otherwise the shape factors would strongly suppress the neutrino 
wave function; thus, \L — Vgt\ < cr^. The physical meaning of the two terms in eq. (38) is 
then clear: A^^.^ is the oscillation phase acquired by the neutrino state over the distance 
L for a pointlike neutrino, whereas A0' takes into account that the wave packet actually 
has a finite size and is the additional phase variation along the wave packet. From eq. (37) 
it is clear that the term A0' can be neglected when the spatial length of the wave packet 
ax is small compared Vg/AE. This coincides with the interaction coherence condition (33) 
considered in Sec. 3.2. We will discuss this condition and the coherence properties of the 
neutrino production and detection processes in more detail in Sec. 5.3. 

Let us now show explicitly that under very general assumptions the neutrino oscillation 
probability takes its standard form (6). We start with the expression for A0 in eq. (37). 
Substituting it into (23), we find 

/,fc(L) = e-^^^ / dtCiL - vg,t) GUL - Vgut) g^^^^-^^-''^*) . (39) 

J — oo 

Let us first neglect the difference between the group velocities of the wave packets describing 
different mass eigenstates, i.e. take Vgi = Vgk = Vg, which also implies = Gi (see Appendix 
A). In this approximation, which neglects the decoherence effects due to the wave packet 
separation, the integral in (39) does not depend on L; this can be readily shown by changing 
the integration variable according to t — ^ (L — Vgt). The same conclusion has already been 
reached in Sec. 3.2 using the momentum-integral representation of Iik{L). The integral in 
eq. (39) is then just the Fourier transform of the squared modulus of the shape factor: 

-/ rfx'|G,(x')pe'^^''*"'. (40) 

^-oo 

It is essentially a coherence factor, which takes into account the effects of suppression of 
the oscillations in the case when the energy and/or momentum uncertainties at neutrino 
production or detection are small enough to allow the determination of the neutrino mass 
(see Sec. 5.3 for a more detailed discussion). If the interaction coherence condition (33) is 
fulfilled, one can replace the oscillating phase factor in the integrand of eq. (40) by unity; 
the resulting integral is then simply equal to 1, as follows from the normalization condition 
(21). Eq. (22) then immediately leads to the standard oscillation probability (6). If, on the 
contrary, the condition opposite to the interaction coherence condition (33) is satisfied, the 
integral (40) is strongly suppressed due to the fast oscillations of the factor exp{iAEikx' /vg) 
in the integrand, leading to the suppression of the oscillations. In the borderline case 
AEik(Tx/vg ~ 1, a partial decoherence occurs. 
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If one now allows for Vgi 7^ Vgk-, then, as direct inspection of the arguments of the 
functions Gi^k in the integrand in eq. (39) shows, the dependence of the integral in this 
equation on L is still negligible if L^Vg/vg <C ax, i.e. if the condition (32) is satisfied. As 
has already been pointed out, this is the condition of the absence of decoherence due to 
the wave packet separation: the distance traveled by neutrinos should be smaller than the 
distance over which the wave packets corresponding to different mass eigenstates separate, 
due to the difference of their group velocities, to such an extent that their effects can no 
longer interfere in the detector.^ If the condition opposite to that in eq. (32) is satisfied, the 
integral in eq. (39) is strongly suppressed because of the lack of overlap between the factors 
Gi{L — Vgit) and Gl.{L — Vg^t) in the integrand (if Iik{L) is written as a momentum-space 
integral, the suppression is due to the fast oscillations of the integrand, see eq. (31)). 

From the above consideration it follows that the factor Iik{L) in the expression for 

the oscillation probability (22) yields the standard oscillation phase factor exp(— z^^^L) 
multiplied by the integral which accounts for possible suppression of the oscillating terms 
due to decoherence caused by the wave packet separation and/or lack of coherence at 
neutrino production or detection. Note that both decoherence mechanisms usually lead 
to exponential suppression of the interference terms in the oscillation probabilities since 
they come from the infinite-limits integrals of fast oscillating functions. The exact form 
of these suppression factors depends on the shape of the wave packets, i.e. is model 
dependent; in particular, for Gaussian and Lorentzian wave packets, these factors are 
~ exp[-(L//coh)^] exp[-(cra;AE)V2t'|] and exp^-L/l^oh) exp{-axAE/2vg), respectively. 

Thus we conclude that the standard oscillation probability is obtained if neutrinos are 
relativistic or quasi- degenerate in mass and the decoherence effects due to the wave packet 
separation or lack of coherence at neutrino production or detection are negligible. No 
unjustified ''same energy" or "same momentum" assumptions are necessary to arrive at this 
result. 

One may wonder why these assumptions are actually so popular in the literature and 
even made their way to some textbooks, though in general there is no good reason to believe 
that the different mass eigenstates have either same momentum or same energy. One possible 
reason could be the simplicity of the derivation of the formula for the oscillation probability 
under these assumptions. However, we believe that simplicity is no justification for using a 
wrong argument to arrive at the correct result. 

^We reiterate that ax is an effective spatial length of the wave packet, which depends on the lengths of 
both the production and detection wave packets axs and axD- The condition (32) therefore automatically 
takes into account possible restoration of coherence at detection, as discussed in Sec. 5.4. 
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4 Another standpoint 



In this section we outline a more general approach to the calculation of the oscillation 
probability, which gives an additional insight into the issues discussed in the present paper. 
Here we just illustrate some points relevant to our discussion rather than presenting a 
complete formalism, which is beyond the scope of our paper. 

The wave functions of the production and detection states \l/f (x, t) and — L) can 

be represented as the integrals over momenta according to (the 1-dimensional versions of) 
eqs. (9) and (15). Inserting these expressions into eq. (18) and performing the integral over 
the coordinate, we obtain 

AabiL, t) = J2 U:,Uu I dq, /z,(g,)e*'''''-*'''(''^)* , (41) 

where 

h^iq)^fnQ-P^)f^'iQ-p'i)■ (42) 

Here the momenta Qi are the integration variables, whereas Pi and p- are, as before, the 
mean momenta of the corresponding wave packets. 

The amplitude (41) is the sum of plane waves corresponding to different momenta and 
different masses. The integration over momenta can be formally substituted by a summation 
to make this point clearer. The oscillation probability is then 



iqiL-iEi{qi)t 



(43) 



The waves with all momenta and masses should be summed up; the resulting expression for 
the oscillation probability includes the interference of these waves. 

The standard approach to the calculation of the oscillation amplitude (41) (or probability 
(43)) is to sum up first the waves with different momenta but the same mass, and then sum 
over the mass eigenstates. In this way first the wave packets corresponding to different mass 
eigenstates are formed, and then the interference of these wave packets is considered. Since 

j dq, /ii(g,)e^''»^-*^»(''')* = Gi{L - Vg.^)e'^^^-'^^^^^^' , (44) 

eq. (41) then directly reproduces the results of Sec. 3. 

Another possibility is to sum up first the plane waves with different masses but equal (or 
related) momenta and then perform the integration over the momenta. In particular, one can 
select the waves with equal energies. Clearly, the final result should not depend on the order 
of summation if no approximations are made, and should be almost independent of this order 
if the approximations are well justified. However, different summation conventions allow 
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different physical interpretations of the result. In what follows we will perform computations 
using the "equal momenta" and "equal energy" summation rules and identify the conditions 
under which they lead to the standard result for the oscillation probability. 

Let us first consider the "equal momenta" summation. Setting qi = p for all i, one can 
write the amplitude (41) as 

Aabim) = J rfp^f/:,[/M/i.(p)e^P^-^^'(^')*. (45) 



From eq. (Al) and the definition (42) of hi{p) it follows that hi{p) = h{p, Ei{p)). We can 

2 . 

i3- 



use this to expand hi in a power series in Amfg • 



h,{p) = h{p, Ei) = h{p, Es) + Ei 



mi=m3 o 



(46) 



dE, 

Inserting this expression into (45), we obtain 

Aat{L,t) = I rfp/.(p,E3(p))e^^^-*^3(^')*5^f/:,f/Me-*[^'(^)-^^(^)l* + ^f,(L,t), (47) 



where 



{p)-E3(p)]t 



Eidhi 



dE, 

It is easy to see that the term ^f^(L,t) is typically very small: 



"ii=m,3 



2Ei{py 



(4J 



Ai{L,t)r.^Au (49) 

where Ai is the first term on the right hand side of eq. (47). Indeed, if the width of the 
effective momentum distribution function h{p) is dp, one has dh/dp ~ h/ap, so that 

Edh Edh Eh Eh 



dE Vg dp VgCp aE 
which immediately leads to (49). Thus, if 



(50) 



'^i^>^, (51) 

^f^(L,t) can safely be neglected. 

Consider now the first term in eq. (47). For a fixed momentum the phase difference 
is A</)j3 = [Ei — E^)t Amf^t/2E-^.^ If the width of the effective momentum distribution 



^This approximation breaks down at very small momenta. Note, however, that the small-p contribution to 
the integral in (47) is strongly suppressed because of the effective momentum distribution function h{p, E3), 
which is strongly peaked at a relativistic momentum p = ps- This justifies using the approximation for A(j)i3 
in (47). 
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function h{p, E^) is small enough, so that the change of the phase within the wave packet 
is small, we can pull the oscillatory factor out of the integral at some effective momentum 
(corresponding to an energy E): 



Aab{L,t) 



tie 



dph{p,E3{p))e'P^-'^'^P^^ 



bie 



G2,{L - Vgst)e 



ip-iL-iEi{p^)t 



(52) 



Here the factor in the square brackets gives the standard oscillation amplitude in the "evo- 
lution in time" approach. Integrating the squared modulus of the amplitude (52) over time 
and using once again eq. (51), one arrives at the standard expression for the oscillation 
probability. 

Thus, we obtain the standard oscillation formula by first summing up the waves with 
equal momenta and different masses and then integrating over the momenta provided that 
the following two conditions are satisfied: 

(i) The variation of the oscillation phase within the wave packet due to the energy spread 
is small: as <^ (27r-E'^/Am|^)L~^; this condition allows one to pull the oscillatory 
factor out of the integral over the momenta, as discussed above. Note that it is 
actually equivalent to the condition of no wave packet separation, eq. (32) (recall that 
Aff, ~ Am^/2_E^ and — Vg/aE)- 

(ii) The momentum distribution functions hi{p) are not too narrow: (Je ^ Am^/2£'. This 
condition, in particular, allows one to neglect the contribution ^^(L,t) in eq. (47). It 
actually ensures that the neutrino wave packet length (Jx — VgjoE is small compared 
to the neutrino oscillation length and is related to the interaction coherence condition, 
as we shall show in Sec. 5.3. 



These conditions for obtaining the standard oscillation formula essentially coincide with the 
conditions found in a different framework in Sec. 3.3. 

To describe the possible decoherence effects due to the separation of the wave packets and 
lack of coherence at neutrino production or detection explicitly, one should lift the conditions 
(i) and (ii) and consider the corresponding corrections to the oscillation amplitude. This is 
discussed in detail in Sec. 5. 

Similarly, we can consider the summation of waves with equal energies and different 
masses, with the subsequent integration over energies (or momenta). Requiring Ei{qi) = 
E^ip) yields qi = ±^/jP~+~Km^^ {i = 1,2), = p. Taking into account that dqi = 
±dp/ ^/l + Am|-/p2, we obtain 



AabiL,t) = J dpY^UliUu 



hi{qi{p),E3{p)) 
y/l + Aml/p^' 



,iqi{p)L-iE3{p)t 



(53) 
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Here it is assumed that 1^3 is the heaviest mass eigenstate, so that Am|j > and no 
singularities appear in the integrand. Note that our change of the integration variables 
Qi ^ p excludes, for i = 1,2, the small regions of momenta qi around zero; this, however, 
introduces only a tiny error, because the main contributions to the integral come from the 
regions around the points qi = Pi, where the functions hi{qi) are strongly peaked. 

Expanding hi{qi{p) , E3{p)) as 

-'h{p,E^{p))+{h{p,E.i{p))-p^i\^^\—f- (54) 



v/1 + Am^,/p2 V dp ) 2p2 

and inserting this expression into (53), we obtain 

Aa,{L, t) = j dp h{p, E3(p))e'^^-^^«(^)* f^aW^e^l^'^^)-^!^ + (55) 



where 



pdh\ Am 



3i 



dp J 2p2 



(56) 



Just as in the previous case, one can show that the amplitude ^f^(L, t) can be neglected if 
the condition (51) is satisfied. Assuming this to be the case and that the variation of the 
oscillation phase within the wave packet due to the momentum spread is small, and taking 
into account that (g^ —p)L = (g^ — ^3)^ ^ Am\^L/2p, one finds 



Aab{L,t) 



2p 



GsiL - Vg3t)e'P'^-'^'^P^^' , (57) 



where p is the average neutrino momentum. This is the standard oscillation amplitude 
multiplied by the effective shape factor of the wave packet of (note that in our current 
approximation we actually neglect the difference between the wave packets of different mass 
eigenstates). Integrating the squared modulus of the amplitude in eq. (57) over time and 
using the normalization condition (21), we again arrive at the standard expression for the 
oscillation probability, just as in the previous case when we first summed the terms with 
equal momenta and different masses and then integrated over the momenta. In deriving 
this result we once again used the conditions (i) and (ii) discussed above. 

Our discussion of the new summation rules for calculating the oscillation probability 
presented here leads to an alternative explanation of why the "same energy" and "same 
momentum" assumptions eventually lead to the correct physical observables, as discussed 
in Sec. 6. 
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5 Quantum- mechanical uncertainty relations and neu- 
trino oscillations 



Neutrino oscillations, being a quantum-mechanical interference phenomenon, owe their very 
existence to quantum-mechanical uncertainty relations. The coordinate-momentum and 
time-energy uncertainty relations are implicated in the oscillations phenomenon in a number 
of ways. First, it is the energy and momentum uncertainties of the emitted neutrino state 
that allow it to be a coherent superposition of the states of well-defined and different mass. 
The same applies to the detection process - for neutrino detection to be coherent, the energy 
and momentum uncertainties inherent in the detection process should be large enough to 
prevent a determination of the absorbed neutrino's mass in this process. The uncertainty 
relations also determine the size of the neutrino wave packets and therefore are crucial to 
the issue of the loss of coherence due to the wave packet separation. In addition, these 
relations are important for understanding how the produced and detected neutrino states 
are disentangled from the accompanying particles. Let us now discuss these issues in more 
detail. We start with 

5.1 Uncertainty relations and disentanglement of neutrino states 

In the majority of analyses of elementary particle processes it is assumed that the energies 
and momenta of all the involved particles have well defined (sharp) values and obey the 
exact conservation laws. However, for this description to be exact, the considered processes 
(and the particles involved) should be completely delocalized in space and in time, whereas 
in reality these processes occur in finite and relatively small spatial volumes and during finite 
time intervals. For this reason, the energy and momenta of all the participating particles 
have intrinsic quantum mechanical uncertainties, and the particles should be described by 
wave packets rather than states of definite momentum - plane waves (see, e.g., [19]). The 
conservation of energy and momentum for these particles is also fulfilled up to these small 
uncertainties. 

This does not, of course, mean that the energy- momentum conservation, which is a 
fundamental law of nature, is violated: it is satisfied exactly when one applies it to all par- 
ticles in the system, including those whose interactions with the particles directly involved 
in the process localize the latter in a given space-time region. Schematically speaking, if we 
consider the process as occurring in a box, the interactions with the walls of the box and the 
contributions of these walls to the energy-momentum balance have to be taken into account. 
In practice, this is never done; however, the resulting inaccuracy of the energy and momen- 
tum conservation as well as the intrinsic quantum-mechanical uncertainties of the energies 
and momenta of the involved particles are usually completely negligible compared to their 
energies and momenta themselves, and therefore can be safely ignored in most processes. 

This is, however, not justified when neutrino oscillations are considered, since the neu- 
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trino energy and momentum uncertainties, as tiny as they are, are crucially important for 
the oscillation phenomenon. In this respect, we believe that the attempts to use the ex- 
act energy-momentum conservation in the analyses of neutrino oscillations are inconsistent. 
In some analyses the exact energy-momentum conservation is assumed for the neutrino 
production and detection processes in order to describe neutrinos as being entangled with 
accompanying particles. The subsequent disentanglement, which is necessary for neutrino 
oscillations to occur, is assumed to be due to the interaction of these accompanying particles 
(such as e.g. electrons or muons produced in decays of charged pions) with medium. This 
localizes those particles and creates the necessary energy and momentum uncertainties for 
the neutrino state. The described approach misses the fact that the parent particles are 
already localized in the neutrino production and detection processes, and so no additional 
disentanglement through the interaction of the accompanying particles with medium is nec- 
essary. Indeed, it is clear that neutrinos produced, for example, in vr^ decays oscillate even 
if the accompanying charged leptons do not interact with medium, i.e. are not "measured". 
The measurement of the flavour of these charged leptons that discriminates between and 
fi^ and makes neutrino oscillations possible is actually provided by the decoherence of the 
charged leptons due to their very large mass difference [7]. 

5.2 What determines the size of the wave packet? 

According to the quantum-mechanical uncertainty relations, the energy and momentum 
uncertainties of a neutrino produced in some process are determined by, correspondingly, 
the time scale of the process and spatial localization of the emitter. These two quantities 
are in general independent; on the other hand, for a free on-shell particle of definite mass 
the dispersion relation E'^ = p'^ + m? immediately leads to 



Since this relation is satisfied for each mass-eigenstate component of the emitted flavour 
state, it must also be satisfied for the state as a whole (provided that the energies and 
momenta of different components as well as their uncertainties are nearly the same, which 
is the case for relativistic or quasi-degenerate neutrinos). Thus, we have an apparently 
paradoxical situation: on the one hand, (Te and cTp should be essentially independent, while 
on the other hand they must satisfy eq. (58). 

The resolution of this paradox comes from the observation that at the time of their 
production neutrinos are actually not on the mass shell and therefore do not satisfy the 
standard dispersion relation. Therefore their energy and momentum uncertainties need 
not satisfy (58). However, as soon as neutrinos move away from their production point 
and propagate distances x such that px ^ 1, they actually go on the mass shell, and 
their energy and momentum uncertainties start obeying eq. (58). This happens because 
the bigger of the two uncertainties shrinks towards the smaller one, so that eq. (58) gets 
fulfilled. Indeed, when the neutrinos go on the mass shell, the standard relativistic dispersion 



EaE = pcTp . 
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relation connecting the energies and momenta of their mass eigenstate components allows 
to determine the less certain of these two quantities through the more certain one, thus 
reducing the uncertainty of the former. As a result, the two uncertainties get related by 
eq. (58), with the one that was smaller at production retaining its value also on the mass 
shell. Note that for neutrino energies in the MeV range neutrinos go on the mass shell as 
soon as they propagate distances x > 10~^° cm from their birthplace. 

Which of the two uncertainties, ap or aE, is actually the smaller one at production? 
Quite generally, this happens to be the energy uncertainty aE- Indeed, consider, e.g., an 
unstable particle, the decay of which produces a neutrino. In reality, such particles are 
always localized in space, so one can consider them to be confined in a box of a linear size 
Ls- The localizing "box" is actually created by the interactions of the particle in question 
with the surrounding particles. Assume first that the average time interval Ts between two 
subsequent collisions of the decaying particle with the walls of the box (more precisely, the 
interval between its collisions with the surrounding particles) is shorter than its lifetime 
r = r~^. Then the energy width of the state produced in the decay is given by the so- 
called collisional broadening, and is actually ~ ^5^^- This width directly gives the neutrino 
energy uncertainty, i.e. aE — Tg^. On the other hand, the neutrino momentum uncertainty 
is ap ^ Lg^. Since Ts is related to Ls through the velocity of the parent particle v as 
Ts — Ls/v, one finds 

as <ap, (59) 

which is actually a consequence of v < 1. 

Consider now the situation when the lifetime of the parent particle is shorter than the 
interval between two nearest collisions with the walls of the box. In this case the decaying 
particle can be considered quasi-free, and the energy uncertainty of the produced neutrino 
is given by the decay width of the parent particle: erg ^ VJ The momentum uncertainty 
of neutrino is then the reciprocal of its coordinate uncertainty a^, which in turn is just 
the distance traveled by neutrino during the decay process: a^ — {p/E)t = {p/E)r~^ ^ 
{p/E)a^^. Thus we find pap ~ EaE, i-e. eq. (58) is approximately satisfied in this case. 
Once again the condition (59) is fulfilled. It can be shown that this inequality is also satisfied 
when neutrinos are produced in collisions rather than in decays of unstable particles [36]. 
Similar arguments apply to the neutrino detection process. Actually, by a^; and ap in the 
above discussion (as well as in the discussion in Sec. 5.3) one should understand the effective 
energy and momentum uncertainties, which depend on the corresponding uncertainties both 
at neutrino production and detection and are dominated, both for the neutrino energy and 
momentum, by the smallest between the production and detection uncertainties. 

Thus, we conclude that the energy uncertainties at neutrino production and detection 
are always smaller than the corresponding momentum uncertainties. This has important 

^This only holds for slow parent particles. In the case of relativistic decaying particles, aE depends on 
the angle between the momenta of the parent particle and of neutrino, see [30] and Sec. 8 below. Still, the 
condition (59) holds in that case as well. 
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implications for the effective neutrino wave packets Gi{L — Vgit). As discussed in Sec. 3.1, 
they are characterized by the effective momentum uncertainty 0"^^ ~ ap = min{aps, CTpo}, 
and similarly for the effective energy uncertainty. Since the wave packets describe prop- 
agating on-shell particles, from the above discussion it follows that the effective energy 
uncertainty characterizing Gi{L — Vgit) coincides with the energy uncertainty at pro- 
duction/detection discussed above: = aE- At the same time, the effective momentum 
uncertainty of the wave packet a'^ has in general nothing to do with the momentum uncer- 
tainties inherent in the neutrino production and detection processes; rather, it is related to 
the energy uncertainty aE through the on-shell condition (58): a'^ = aE/vg- Because the 
spatial length of the wave packets is a^ ~ V'^p^? ^^is in turn means that 

CTx^ Vg/aE- (60) 

Thus, we conclude that the spatial length ax of the wave packets describing the propa- 
gating neutrino states is always determined by the energy uncertainty at neutrino produc- 
tion/detection as the smaller one between cXp and aE- This is in accord with the known fact 
that for stationary neutrino sources (for which = 0) the neutrino coherence length is 
infinite [35, 37]. On the other hand, as will be shown in the next subsection, the localization 
conditions for the neutrino production and detection processes, which have to be fulfilled 
for these processes to be coherent, are always determined by the corresponding momentum 
uncertainties. 

An interesting case in which the on-shell condition dp ^ 0"e is strongly violated both 
at production and detection is the proposed Mossbauer neutrino experiment, for which 
aE ~ lO""*^"^ eV and ap ~ 10 keV are expected [1, (3], so that aE ~ 10~^^ap. We discuss some 
important implications of this large disparity between aE and ap in Sec. 5.4. 

Our final comment in this section is on the case when ap ~ aE/vg at neutrino production 
and detection, so that ax ~ a'^. While the localization condition requires relatively large ap 
{ap ^ Am^/2p) for the emitted and detected neutrino states to be coherent superpositions 
of mass eigenstates (see Sec. 5.3), the condition of no decoherence due to the wave packet 
separation, on the contrary, requires long wave packets, i.e. relatively small ap. Is there any 
clash between these two requirements? By combining the two conditions we find Am^/2p ^ 
ap <S {vg/ Avg)L~^, which can only be satisfied if 

AmV2j9 < {vg/Avg)L-^ . (61) 

This can be rewritten as the following condition on the baseline L: 

27r-^«-^. (62) 

^osc 9 

Since Vg/Avg ^ 1, this condition is expected to be satisfied with a large margin in any 
experiment which intends to detect neutrino oscillations: if it were violated, neutrino oscil- 
lations would have been averaged out because of the very large oscillation phase (except for 
unrealistically good experimental energy resolution 6E/E < Avg/vg ~ Am^/2£'^). 
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5.3 Coherence of the produced and detected neutrino states 

In order for a neutrino state produced in a charged- current weak interaction process to be 
a coherent superposition of different neutrino mass eigenstates, it should be in principle 
impossible to determine which mass eigenstate has been emitted. This means that the in- 
trinsic quantum-mechanical uncertainty of the squared mass of the emitted neutrino state 
(7^2 must be larger than the difference Am^ of the squared masses of different neutrino mass 
eigenstates [21, 23]: (7^2 ^ Am^. Conversely, if am? -C Am^, one can determine which mass 
eigenstate has been emitted, i.e. the coherence of different mass eigenstates is destroyed. 
This situation is quite similar to that with the electron interference in double slit experi- 
ments: If there is no way to find out which slit the detected electron has passed through, 
the detection probability will exhibit an interference pattern, but if such a determination is 
possible, the interference pattern will be washed out. 

Assume that by measuring the energies and momenta of the other particles involved 
in the production (or detection) process we can determine the energy E and momentum 
p of the emitted or absorbed neutrino state, and that the intrinsic quantum-mechanical 
uncertainties of these quantities are aE and CTp. From the energy - momentum relation 
E"^ = -^m? we can then infer the squared mass of the neutrino state with the uncertainty 
am? = [(2-E(T£;)^+(2p(Tp)^]^''^, where it is assumed that aE and ap are uncorrelated. Therefore 
the condition that the neutrino state be emitted or absorbed as a coherent superposition of 
different mass eigenstates is [21, 23] 

= [{2EaEf + (2|)crp)^] ^'^ > Am^ . (63) 

This condition has a simple physical meaning. As follows from eq. (59), for relativistic 
neutrinos the second term in the square brackets in eq. (63) exceeds the first one, and 
therefore this condition essentially reduces to 2pap ^ Am^, or 

cTp > AmV2j9 = 27r//osc • (64) 

This is the so-called localization condition for neutrino production and detection. Indeed, 
since the effective production/detection momentum uncertainty ap is related to the coordi- 
nate uncertainties of the neutrino emitter and absorber o"^'^ and a^'^ by 
ap ~ mm{(cj^°^)-\ (a^')~^}, eq. (64) yields 

a^°^«/osc, ^^^«/osc. (65) 

This is nothing but the obvious requirement that the neutrino production and detection 
processes be localized in spatial regions that are small compared to the neutrino oscillation 
length; if it is violated, neutrino oscillations will be averaged out upon the integration over 
the neutrino production and detection coordinates in, respectively, the neutrino emitter and 
absorber. Such a washout of the oscillations is equivalent to decoherence. 

In real situations one always deals with large ensembles of neutrino emitters, and the 
detectors also consist of a large number of particles. Therefore, in calculating the observable 
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quantities - the numbers of the neutrino detection events - one always has to integrate over 
the macroscopic volumes of the neutrino source and detector. In some situations (e.g., for 
solar or reactor neutrinos) the source and detector are much larger than the localization 
domains of the wave functions of individual neutrino emitters and absorbers. In these cases 
the integration over the source and detector volumes modifies the localization conditions: 
instead of depending on the spatial sizes of the individual neutrino emitter and absorber, a^'^ 
and (T^^, they contain the macroscopic lengths of the source and detector in the direction of 
the neutrino beam, Ls, and L^. In other words, a necessary condition for the observability 
of neutrino oscillations is 

Ls,D < /osc , (66) 

which can be much more restrictive than the conditions in eq. (65). If this condition is 
violated, neutrino oscillations are averaged out. 

Let us now discuss the coherence condition for the neutrino production and detection 
processes from a slightly different perspective. Quite generally, this condition can be formu- 
lated as a requirement that the energy and momentum uncertainties inherent in the neutrino 
production and detection processes be much larger than, correspondingly, the energy and 
momentum differences of different mass eigenstates AE and Ap: 

\AE\<^aE, \Ap\<^(Tp. (67) 

As mentioned above, aE is the effective energy uncertainty, which depends on the energy 
uncertainties both at production and detection and is dominated by the smallest between 
them, and similarly for the momentum uncertainty ap. It is now easy to see that if both 
the conditions in eq. (67) are fulfilled, eq. (63) is satisfied as well. 

Let us now find the relationship between the conditions (67) and the one in eq. (33), 
which we called the interaction coherence condition in Sec. 3.2. As was pointed out in the 
previous subsection, the length of the wave packet is determined by the effective energy 
uncertainty aE- — Vg/aE- Therefore one can rewrite the condition (33) as 



\AE\ 
(Te 



< 1 . (6^ 



This is nothing but the first condition in (67). To study the condition (33) further, let us 
rewrite (68), by making use of eq. (34), as 



VgAp ^ Am^ 



aE 2Ea£ 



< 1 . (69) 



Now, there are essentially two possibilities: 

(A) VgAp is not approximately equal to —Aim?/2E, and 



^Indeed, by making use of the relation Am^ ~ 2EAE — 2pAp (which follows from — — 

1 /2 

and is equivalent to (34)), one can rewrite (63) as [{2EaE)^ + {2pap)^] > \2EAE - 2pAp\. From this 
representation of (63) the above statement is evident. 
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(B) VgAp ^ -Am'^/2E. 

Let us now consider these cases in turn. 

(A) . In this case there are no cancellations between the two terms on the l.h.s. of 
eq. (69). From kinematics considerations it follows that case A is realized when the energy 
release at neutrino production is comparable to (or much larger than) the masses of the 
accompanying particles, as e.g. in the decays tt"*" — * /i+ + z/^j and tt"*" ^ e"*" + z/g- In this case 
AE ~ Ap. In the absence of cancellations between the two terms on the l.h.s. of (69) the 
condition (69) implies 

« 1 , « 1 • (70) 

(Te 2EaE 

Since o"p > (Je-, from the first strong inequality in (70) we find that for relativistic neutrinos 
\Ap\/ap <ti 1, which is the second condition in (67). Thus, in case A the condition (33) 
enforces those in eq. (67). 

Because aE — Vg/a^ and /qsc = 47rp/Am^, the second condition in (70) is equivalent to 
(^x ^osc; i-e. in this case the length of the neutrino wave packet is small compared to the 
oscillation length. 

(B) . Case B corresponds to a strong cancellation between the two terms on the l.h.s. 
of eq. (69) (and between the two terms on the r.h.s. of eq. (34)). From eq. (34) it follows 
that in this case \AE\ -C |Ap|, which is close to the "same energy" situation. Case B is 
realized when the energy release in the neutrino production process is small compared to 
the masses of some of the accompanying particles; examples are nuclear (3 decay and the 
decay vr"*" — > vr'' + e"*" + i/g- Case B can be further subdivided into two cases: 

• Case Bi: both terms on the l.h.s. of (69) are small, i.e. the conditions (70) are satisfied. 
The situation with coherence of neutrino production and detection in this case is the 
same as in case A, and therefore case Bi does not require a separate analysis. 

• Case B2: both terms on the l.h.s. of (69) are not small, i.e. 

> 1 ^ > 1 (71) 

and only their sum is small. Note that the signs > here can even be replaced by ^ 
provided that the cancellation between the two terms in (69) is almost exact. Since 
ce < CTp (and may also be <C ap), the first condition in (71) does not necessarily mean 
that the momentum differences of different mass eigenstates exceed the momentum 
uncertainties at production or detection, i.e. in general it does not contradict the 
second condition in (67); on the other hand, the condition \Ap\ <^ ap is not in this 
case automatically enforced by (33). 

The second inequality in (71) implies that the length of the neutrino wave packet 
is not small compared to the oscillation length. This is the case, in particular, for 
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the proposed Mossbauer neutrino experiments [4, 6]. Note that the fact that the 
length of the wave packet can in some cases be larger (or even much larger) than 
the oscillation length does not mean that the oscillations cannot be observed in those 
cases. Indeed, from eq. (37) it follows that, for a fixed L, the variation of the phase 
difference Ac/) due to the variation of the quantity \L — Vgt\ between and ax does 
not exceed {AE/vg)ax. The condition (68) (which is equivalent to the one in eq. (33)) 
ensures that this variation is negligible. This situation actually corresponds to a quasi- 
stationary case, when the oscillatory pattern depends on the spatial coordinate but 
practically does not change with time. Therefore, if the sizes of the neutrino source 
and detector are small compared to the oscillation length, non-averaged oscillations 
will be observed despite the very large length of the wave packet. 

Thus, we conclude that, in the wave packet approach, in cases A and Bi the condition 
(33) guarantees the coherence of the neutrino production and detection processes. At the 
same time, in case B2 characterized by eq. (71), the condition (33) is not sufficient to ensure 
the coherence of neutrino production and detection and one has to supplement it by the 
second condition in (67), which in this case does not automatically follow from the wave 
packet formalism. Note that a strong cancellation between the two terms in (69) implies 
|Ap| ^ Am^/2p, and therefore the second condition in (67) coincides in this case with the 
localization condition (64)). 

The above observations underline an important difference between the the wave packet 
approach and the more consistent quantum field theoretical (QFT) treatment of neutrino 
oscillations, in which neutrino production, propagation and detection are considered as a 
single process (see, e.g., [38, 39, 40, 23, 37, 41, 42, 26, 27, 36, 28, 43, 44, 4]). The wave 
packet approach does not take into account the neutrino production and detection processes, 
except by assigning to the neutrino state a momentum uncertainty, which is supposed to 
be determined by these processes. In particular, the wave packet picture assumes the mass- 
eigenstate components of the fiavour neutrino states to be always on the mass shell, so 
that their energy and momentum uncertainties are always related by eq. (58). As discussed 
above, the wave packet picture is quite adequate in cases A and Bi, but it needs to be 
supplemented by the information on the source and detector localization in case B2. At the 
same time, in the QFT-based treatment of neutrino oscillations, a factor that takes into 
account the neutrino emitter and absorber localization properties emerges automatically in 
the expression for the oscillation probability in all cases. 

This being said, we should add that, although case B discussed above is by no means 
rare, we are currently aware of only one instance when the situation in that case corre- 
sponds to the subcase B2, i.e. is characterized by the conditions (71): this is the proposed 
Mossbauer neutrino experiment. To the best of our knowledge, for neutrinos from con- 
ventional sources the fulfillment of (33) always guarantees the coherence of the neutrino 
production and detection processes. It is for this reason that we dubbed it the interaction 
coherence condition. 
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5.4 Wave packet separation and restoration of coherence at de- 
tection 

Let us now assume that a neutrino flavour state was produced coherently in a weak inter- 
action process and consider its propagation. The wave packet describing a flavour state 
is a superposition of the wave packets corresponding to different mass eigenstates. Since 
the latter propagate with different group velocities, after some time tcoh they will separate 
in space and will no longer overlap. If the spatial length of the wave packet is ax-, this 
time is tcoh — o'x/ ^Vg. The distance /coh that the neutrino state travels during this time is 
/coh — Vg{ax/Avg). If the distance L between the neutrino emission and detection points 
is small compared to the coherence length, i.e. if the condition (32) is satisfied, then the 
coherence of the neutrino state is preserved, and neutrino oscillations can be observed. If, 
on the contrary, L ~ l^oh or L S> Icoh, partial or full decoherence should take place. 

If the wave packet length in the above discussion is assumed to be fully determined 
by the emission process, then this is not the full story yet: even if the coherence is lost 
on the way from the source to the detector, it still may be restored in the detector if the 
detection process is characterized by high enough energy resolution [35]. According to the 
quantum-mechanical time-energy uncertainty relation, high energy resolution requires the 
detection process to last sufficiently long; in that case, the different wave packets may arrive 
during the detection time interval and interfere in the detector (or, more precisely, the effects 
produced in the detector by different wave packets interfere). One can take into account 
this possible restoration of coherence in the detector by considering to be an effective 
length of the wave packet, which exceeds the length of the emitted wave packet and takes the 
coherence of the detection process into account. This effective wave packet length is actually 
the length that characterizes the shape factors Gi{L — Vgit), as discussed in Sees. 3.1 and 
5.2. With ax being the effective wave packet length, the coherence condition (32) includes 
the effects of possible restoration of coherence at detection. 

An interesting example of coherence restoration at detection, which resolves another 
paradox of neutrino oscillations, is the proposed Mossbauer neutrino experiment. In this 
case aE ~ 10~^^ eV and cxp ~ 10 keV are expected [4, 6]. In such an experiment neutrinos 
are produced coherently due to their large momentum uncertainty [1]. However, as soon as 
the emitted neutrino goes on the mass shell, its momentum uncertainty shrinks to satisfy 
eq. (58), i.e. essentially becomes equal to the tiny energy uncertainty. Therefore for on- 
shell Mossbauer neutrinos a^^ ~ 10"^^ eV ^ ^m?/2E ^ 10"'' eV, i.e. the momentum 
uncertainty is much smaller than the difference of the momenta of different mass eigenstates. 
This means that coherence of different mass eigenstates in momentum space is lost. However, 
the fact that in this case both the energy and momentum uncertainties of the propagating 
neutrino state are much smaller than iSvp? jlE does not mean that oscillations cannot be 
observed. In fact, it has been shown in [ I] that the Mossbauer neutrinos should exhibit 
the usual oscillations. The resolution of the paradox lies in the detection process: the large 
momentum uncertainty at detection, (Xp ~ 10 keV ^ Am^/2i?, restores the coherence by 
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allowing the different mass eigenstates composing the ffavour neutrino state to be absorbed 
coherently. 

It is well known that coherence plays a crucial role in observability of neutrino oscilla- 
tions. It is interesting to note, however, that even non-observation of neutrino oscillations 
at baselines that are much shorter than the oscillation length is a consequence of and a firm 
evidence for coherence of the neutrino emission and detection processes: if it were broken 
(i.e. if the different neutrino mass eigenstates were emitted and absorbed incoherently), the 
survival probability of neutrinos of a given flavour, instead of being practically equal to 1, 
would correspond to averaged neutrino oscillations. 

6 When is the stationary source approximation justi- 
fied? 

It has been pointed out in [23, 35] and elaborated and exploited in [14] that for stationary 
neutrino sources the following two situations are physically indistinguishable: 

(a) A beam of plane-wave neutrinos, each with a definite energy E and with an overall 
energy spectrum $(£"); 

(b) A beam of neutrinos represented by wave packets, each of them having the energy 
distribution function f{E) such that \f{E)\'^ = ^{E). 

As was stressed in [14], this actually follows from the fact that in stationary situations 
the spectrum $(£') fully determines the neutrino density matrix and therefore contains the 
complete information on the neutrino system. 

This, in fact, gives an alternative explanation of why the "same energy" approach, 
though based on an incorrect assumption, leads to the correct result. It has been shown 
in Sec. 4 that, within the proper wave packet formalism, one can choose to sum up first 
the states of different mass but the same energy, and then integrate over the energy (or 
momentum) distributions described by the effective energy or momentum shape factors of 
the wave packets h{E) or h{p). In the light of the physical equivalence of the situations (a) 
and (6), it is obvious that the integration over the spectrum of neutrinos, which is inherent 
in any calculation of the event numbers, leads to the same result as the integration over the 
energy spread within the wave packets (provided that the corresponding energy distributions 
coincide). Thus, the "same energy" assumption, though by itself incorrect, leads to the 
correct number of events upon the integration over the neutrino energy spectrum. The same 
is true for the "same momentum" assumption. This actually means that the wave packet 
description becomes unnecessary in stationary situations, when the temporal structure of 
the neutrino emission and detection processes is irrelevant and the complete information on 
neutrinos is contained in their spectrum $(-£"), as was first pointed out in [14]. 
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Let us derive the results of stationary source approximation in terms of the wave packet 
picture described in this paper. We start with eq. (41) for the oscillation amplitude with 
hi{p) defined in (42). Notice that the momentum distribution functions ff{p — Pi) and 
fFiP ~ P'i) introduced in eqs. (9) and (15) do not depend on time, and furthermore they 
are defined for all moments t from — oo to +oo. The only time dependence in (9) is in 
the form of the plane waves in the integrand. This is precisely what corresponds to the 
stationarity condition: the source has no special time feature, and there is no tagging of 
neutrino emission and detection times. 

From eqs. (41) and (42) we find 

\Aatmt)\' = J2U:,U^UakU;, J dp J dp'ff{p-p,)ft{p'-Pk)fF*{p-p'^fk{p'-p'k) 

i,k 

^(,i{p-v')L ^-i[E,(p)-Et,{p')]t ^ ^^2) 

The integration over time is trivial: 

/ + CXD 
dt e-[^«(p)-^^(p')]* = 2tx5[E,{p) - Ek{p')] , (73) 
-oo 

which means that only the waves with equal energies interfere. We stress once again that 
this is a consequence of the fact that no time structure appears in the detection and produc- 
tion processes, which is reflected in the time independence of the momentum distribution 
functions ff{p—pi) and fP{p — p'i)-, and in the integration over the infinite interval of time. 
The equality Ei{p) = Ek{p') leads to 

p-p' = Aml/2p (74) 

to the leading order in the momentum difference. The (5-function (73) can be used to remove 
one of the momentum integrations in (72), so that we finally obtain 

, ik^ 



/H-OO r) p 

dt\Aat{L,t)\' = — dp\f{p-p)\V{p-p')\'J2u:,U,.U:,Uake ^ -V 
-oo J ■ , 

^ t,k 

(75) 

Here we have neglected the dependence of the shape factors on the neutrino mass. Replacing 
the integration over momenta by the integration over energies, we can rewrite the oscillation 
probability as 

/+00 9 /• 

dt\Aab{L,t)\^ = — / dE<l>{E)R{E)P,,{E,L), (76) 
-oo J 

where 

_ _ Am?, 

P,,iE, L) = J2 K^uUt^Ua.e-^^'^ (77) 
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is the standard expression for the oscillation probability, = {/^{E — E)]'^ is the energy 

spectrum of the source, and R{E) = \ f^{E — E')\'^ is the resolution function of the detector 
(note that we have substituted the momentum dependence of these quantities by the energy 
dependence using the standard on-shell dispersion relation). 

The expression in eq. (76) corresponds to the stationary source approximation: the 
oscillation probability is calculated as an incoherent sum of the oscillation probabilities, 
computed for the same-energy plane waves, over all energies. 

Notice that we have performed integration over the spatial coordinate at the level of the 
amplitude and over time at the probability level. Apparently, such an asymmetry of space 
and time integrations is not justified from the QFT point of view. In QFT computations the 
integration over time is performed in the amplitude, and this leads (in the standard setup) 
to the delta function which expresses the conservation of energy in the interaction process. 
To match our picture with that of QFT we need to consider the detection process and take 
into account the energies of all the particles that participate in the process. Suppose that 
the algebraic sum of the energies of all the accompanying particles (taken with the "-" sign 
for all incoming particles the "+" sign for the outgoing ones) is Ed- Then instead of (73) 
we will have in the probability 

P + OO P + OO 



dt e-^[^«(p)-^°l* / dt' e^[^^(f')-^^l*' = 4tt^ S[E,{p) - Ek{p')] S[Ei{p) - Ed] , (78) 



where the second J-function on the right hand side reflects the energy conservation in the 
detection process. Using (78) we again obtain the "same energy" interference, as before. 

In the above calculation we have not introduced any time structure at detection and 
performed the integration over t from — oo to +00. In reality, certain time scales are always 
involved in the detection processes (even if we do not perform any time tagging in the 
emission process). For example, we can measure with some accuracy the appearance time 
of a charged lepton produced by the neutrino capture in the detection process. In this case, 
the neutrino detection state will have a time dependence: 



where has a peak at to of a width at that is determined by the accuracy of the mea- 
surement of the time of neutrino detection. Since in practice the spatial characteristics of 
neutrino detection do not change with time, the dependences of on x and t factorize: 
\l/f = "^xii^ ~ L)^ti{'t " ^o)- Integrating over time in the amplitude, we will have 



dx / dt *(a: -L,t- to) e^P^-^(^'-^°)* = 2'Ke'P^f^*{p - p[) f^^*{Ei - Ed) , (80) 






(79) 




where 
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is the Fourier transform of — to)- As we have mentioned, — to) has a peak 

of the width 0"^ at t = to- Taking for an estimate at ~ 10"^ s (which is probably the 
best currently achievable time resolution), we obtain 5E ~ a^^ ~ 10"^ eV. This is many 
orders of magnitude smaller than the typical energy resolution in the oscillation experiments. 
Therefore one can substitute fu{Ei — Ed) \/2ti 5{Ei — Ed)-, which brings us back to our 
previous consideration. 

7 When can neutrino oscillations be described by pro- 
duction and detection independent probabilities? 

In most analyses of neutrino oscillations it is assumed that the oscillations can be described 
by universal, i.e. production and detection process independent probabilities. In other 
words, it is assumed that by specifying the flavour of the initially produced neutrino state, 
its energy and the distance between the neutrino source and detector, one fully determines 
the probability of finding neutrinos of all flavours at the detector site (for known neutrino 
mass squared differences and leptonic mixing matrix). The standard formula for neutrino 
oscillations in vacuum, eq. (6), is actually based on this assumption. Such an approach is 
very often well justified, but certainly not in all cases. It is, therefore, interesting to study 
the applicability limits and the accuracy of this approximation. 

A natural framework for this is that of QFT, which provides the most consistent ap- 
proach to neutrino oscillations. In this method the neutrino production, propagation and 
detection are considered as a single process with neutrinos in the intermediate state. This 
allows one to avoid any discussion of the properties of the neutrino wave packets since neu- 
trinos are actually described by propagators rather than by wave functions. The properties 
of neutrinos in the intermediate state are fully determined by those of the "external" parti- 
cles, i.e. of all the other particles that are involved in the neutrino production and detection 
processes. The wave functions of these external particles have to be specified. Usually, these 
particles are assumed to be described by wave packets; for this reason the QFT-based treat- 
ment is often called the "external wave packets" approach [32], as opposed to the usual, 
or "internal wave packets" one, which was discussed in Sees. 3, 4 and 6 and which does 
not include neutrino production and detection processes. The results of the QFT-based 
approach turn out to be similar, but not identical, to those of a simple wave packet one; in 
particular, possible violations of the on-shell relation (58) between the neutrino energy and 
momentum uncertainties is now automatically taken into account. Moreover, the values of 
these uncertainties, which specify the properties of the neutrino wave packet in the "internal 
wave packets" approach and which have to be estimated in that approach, are now directly 
derived from the properties of the external particles. 

The results of the QFT approach can be summarized as follows. For neutrinos propagat- 
ing macroscopic distances the overall probability of the production-propagation-detection 
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process for relativistic or quasi-degenerate neutrinos can to a very good accuracy be rep- 
resented as a product of the individual probabilities of neutrino production, propagation 
(including oscillations) and detection.^ The oscillation probability, however, is not in gen- 
eral independent of production and detection processes, which means that the factorizability 
of the probability of the entire production - propagation - detection process and the univer- 
sality of the oscillation probability (or lack thereof) are in general independent issues. The 
oscillation probability can be generically represented as 

P{va - L)=J2 K^uUakU;^ e-^^^ 5coh(^//?f ) 5p/D ( Am^Ja™,. ) . (82) 

i,k 

Here S'p/D(Am^;,/(Tm2) and Scoh{L/l^^) are, respectively, the production/detection and prop- 
agation coherence factors, which account for possible suppression of the oscillations due to 
the lack of coherence at neutrino production or detection and due to the wave packet sep- 
aration. They are both equal to unity at zero argument and quickly decrease (typically 
exponentially) when their arguments become large. The simple "internal wave packets" 
approach leads to an expression for the oscillation probability that is similar in form to that 
in eq. (82) but with the argument of S'p/d replaced by (AEik/cTE). As discussed in Sec. 5.3, 
the two production/detection coherence conditions, Am^^ <^ (7^2 and AEik <C cte, are 
essentially equivalent except when the conditions (71) are satisfied. Note that in the QFT- 
based framework the form of the propagation and production/detection coherence factors 
Scoh{L/li^^) and S'p/d (Am^^/cXma) in eq. (82) is derived from the properties of the external 
particles and of the detection and production processes. 

For the oscillation probability to be independent of the processes of neutrino production 
and detection, the following conditions have to be fulfilled: 

(i) The neutrino emission and absorption should be coherent, and decoherence effects due 
to wave packet separation should be negligible; 

(ii) The energy release in the production and detection reactions should be large compared 
to the neutrino mass (or compared to the mass differences). 

The necessity of (i) is clear from the discussion above: if this condition is fulfilled, both 
the coherence factors in eq. (82) are equal to unity, and the standard neutrino oscillation 
formula is recovered. If, on the contrary, (i) is violated, the oscillations will in general suffer 
from the production and detection dependent decoherence effects (note, however, that if the 
decoherence is complete, the oscillation probability will still have a universal form as it will 
correspond to averaged oscillations in that case). As to the condition (ii), it ensures that 

notable exception from this rule is the case of the Mossbauer effect with neutrinos, in which the 
probabilities of neutrino production and detection do not factorize but are instead entangled with each 
other [4]. Still, even in this case, the oscillation (actually, Ve survival) probability can to a very good 
accuracy be factored out of the expression for the overall probability of the process. 
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the production and detection probabilities are essentially the same for all mass-eigenstate 
components of the emitted or detected flavour neutrino states (modulo the different values of 
|[/aip); if this condition is violated, the phase space available in the production or detection 
process will depend on the mass of the participating neutrino mass eigenstate, and the 
mass-eigenstate composition of the flavour eigenstates will no longer be given by the simple 
formula (1). 

8 Wave packet approach and Lorentz invariance 

Let us start the discussion of the Lorentz invariance issues with considering yet another 
apparent paradox. In Sec. 5.2 it was pointed out that the energy uncertainty of neutrinos 
produced in decays of free or quasi-free unstable particles at rest is essentially given by 
the parent particle's decay width, i.e. aE — L. Correspondingly, the spatial length of the 
neutrino wave packet is 



where r is the mean lifetime of the decaying particle. 

Consider now the situation when in the laboratory frame the parent particle moves with 
a speed u. One might then expect that the neutrino energy uncertainty in the laboratory 
frame a'^ will be given by the decay width of the unstable particle in this frame F' = r/7„, 
where 7„ = 1/ yl — v? is the Lorentz factor of the parent particle. The length of the neutrino 
wave packet a'^ is then expected to be v'g/V ~ l/F' = 7M/r, which means that a'^ should be 
increased by a factor of 7„ compared to a^- On the other hand, the transition from the rest 
frame of the unstable particle to the laboratory frame is described by a Lorentz boost; if 
the parent particle is boosted in the direction of the neutrino momentum, then the length of 
the neutrino wave packet should be reduced due to the Lorentz contraction, not increased! 

The resolution of this paradox comes from the observation that the decay process takes a 
finite time r' = l/F', during which the decaying particle propagates some distance. Consider 
the situation when the neutrino is emitted along the momentum of the parent particle. 
During the decay time the parent particle moves over distance / = ut' in the direction of 
the motion of the neutrino, thus reducing the length of the emitted neutrino wave packet: 



^ Vgl as ^ Vg/V = VgT, 



(83) 




(84) 



where in the last equality we used the relativistic law of addition of velocities 



Vg+U 



(85) 



1 + VgU 



Comparing eq. (84) with eq. (83), we find 



a. 



X 



(86) 



7„(1 + VgU) ' 
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which is just the standard relativistic law of transformation of lengths in the direction of 

^9 — 



the motion of objects. ^'^ For relativistic neutrinos f „ ~ ~ 1 and 



+ u 

7,(l + ygM) ~7„(l + u) = ^/^ , (87) 

— u 

i.e. the neutrino wave packets are contracted when neutrinos are emitted in the direction 
of the motion of the parent particles and dilated when they are emitted in the opposite 
direction {u < 0), as it should be. This result was previously obtained on the basis of 
different arguments in [30]. 

Let us now discuss the behaviour of the oscillation probability with respect to the 
Lorentz boosts. We first analyze the standard expression for the oscillation probability (6) 
(our discussion here mostly follows that in [15]; a more general though less direct proof will 
be given below). Consider a boost with the speed u of the neutrino source - detector system 
in the direction of the neutrino motion (the generalization to the boost in an arbitrary 
direction is straightforward). Under this boost the distance L traveled by the neutrino 
between its emission and absorption points, the propagation time t and the neutrino energy 
and momentum transform according to 

L' = 7„(L + ut) , t' = 7,(t + uL) , (88) 
E' = ^^{E + up), p' = ^u{p + uE). (89) 

Recall now that the standard formula for the oscillation probability results when the pro- 
duction/detection and propagation coherence conditions are satisfied; for neutrinos from 
conventional sources the former condition includes the requirement that the neutrino wave 
packet length be small compared to the oscillation length /qsc (Mossbauer neutrinos, for 
which this condition is not satisfied, represent a special case, to be considered below). Since 
the baselines of interest in neutrino oscillation experiments satisfy L > Zqsc, the condition 
o'x <^ ^osc also means that the length of the neutrino wave packets is negligible compared 
to the baseline, so that one can consider neutrinos pointlike and set L = Vgt. Substituting 
this into the first equality in eq. (88) yields L' = 7mL(1 + u/vg). On the other hand, since 
= p/E, the second equality in eq. (89) can be rewritten as p' = 7„j'(l + u/vg). Thus 

^ = -, (90) 
p p 

which proves Lorentz invariance of the standard oscillation probability (6). 



^°Upon a boost with the speed v from its rest frame to a frame K , the length of an object in the direction 
of the boost cP is contracted by the factor 71, = l/Vl — v^: d = (P/"fv In the reference frame K' which 
moves in the same direction as K with the speed u with respect to the latter, one has d' = d'^/jv', where 
v' = {v + m)/(1 + vu). From the latter equality one finds j^i = 7„7«(1 + uv). This gives d' = d^y/jyi = 
(i/[7„(l + vu)], to be compared with eq. (86). 
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Note that the first equahty in (89) can be written as E' = -E'7u(l + uvg), and therefore 
eq. (86) imphes that the product a^E is invariant under Lorentz transformations. This 
result was previously obtained in [30] from different considerations. 

We shall now present an alternative (and more general) proof of Lorentz invariance 
of the standard oscillation probability. The phase difference A0 defined in eq. (24) is 
Lorentz invariant; hence so is the expression for A0 in eq. (37), which was obtained from 
(24) assuming neutrinos to be relativistic or quasi-degenerate. (Recall that we adopt this 
assumption throughout this paper, and that the standard expression for the oscillation 
probability (6) is only valid under this assumption). The two terms in eq. (37) are not in 
general separately Lorentz invariant, though their sum is. However, in the situations when 
the second term is negligibly small in all Lorentz frames, the first term must be Lorentz 
invariant by itself; the standard oscillation probability (6) is then Lorentz invariant. As 
discussed in Sees. 3.2 and 5.3, the second term in eq. (37) is responsible for possible violation 
of the production/detection coherence condition, and can be discarded when this condition 
is satisfied. From eq. (37) it follows that this term vanishes, in particular, in the limit of 
pointlike neutrinos, L = Vgt, which leads to Lorentz invariance of the standard oscillation 
probability, as discussed above. In addition, this term vanishes in the limit of vanishingly 
small AE even if the condition L = Vgt is not satisfied; this case is realized for Mossbauer 
neutrinos [4, 6]), whose oscillation probability is therefore also Lorentz invariant. We add 
here that the condition L = Vgt is Lorentz invariant, and so is the condition AE ^ 
for Mossbauer neutrinos; the former follows directly from eq. (88), while the latter is a 
consequence of the fact that Mossbauer neutrinos are emitted and absorbed recoillessly in 
any Lorentz frame. 

How about the general case when the production/detection and propagation coherence 
conditions are not necessarily satisfied? The oscillation probability must be Lorentz invari- 
ant in that case as well. In the QFT-based approach, this comes out automatically since the 
calculations are based on Lorentz covariant Feynman rules. In contrast to this, calculations 
in the wave packet approach are not manifestly Lorentz covariant; however, as we shall show 
now, they still lead to the Lorentz invariant expression for the oscillation probability. 

Let us first demonstrate that the integral over time of the product of any two wave 
packet functions, Gi{L — Vgit) and Gl.{L — Vgkt), is Lorentz invariant. We write the effective 
wave packets as Gi = |Gj|e*^'. From the normalization condition (21) it follows that the 
quantities \Gi\'^dt are Lorentz invariant for all i, i.e. 

\Gi\^dt = inv , \G k\'^ dt = inv . (91) 

Taking the square roots of these two equalities and multiplying the results, we find 

\Gi\\Gk\dt = inv. (92) 

Next, we consider the complex phase factors e^^'. They are Lorentz invariant because so 
must be the phases pi. The latter is a consequence of the fact that for any physical quantity 



35 



that can be expanded in a nontrivial power series the expansion parameter must be Lorentz 
invariant - otherwise different terms in the expansion would transform differently under the 
Lorentz transformations. From the invariance of e*''' and eq. (92) it follows that the integral 
over time of the product of any two wave packet functions, Gi{L — Vgit) and G*f,{L — Vgkt), 
is Lorentz invariant, as advertised. Since so is the phase difference A0, the quantities 
Iik{L) defined in eq. (23) are Lorentz invariant as well. Lorentz invariance of the oscillation 
probabilities Pab{L) then follows immediately from eq. (22). 

Before closing the discussion of the Lorentz invariance issues, let us note an interesting 
feature of the wave packets Gi{L — Vgit). As discussed in Sec. 3.1, their moduli decrease 
rapidly when \L — Vgit\ exceeds the effective length of the wave packet a^- Therefore \Gi\ 
should depend on (L — Vgit) through the ratio 

Ri = . (93) 



Under the Lorentz transformations we have 



L' - v'g,t' = 7, 



Vgi + u 



(L + ut) - (t + uL) 



7„(1 + VgiU) 



""''^ (94) 



From this formula and the properties of ax under Lorentz transformations (86) it follows 
that the ratio (93) is Lorentz invariant. Note that, in addition to their dependence on R^, 
\Gi{L — Vgit)\ depend on a^i and Vgi through their normalization factors. 



9 Remarks on the wave packet approach 

We add here a few comments on the wave packet approach to neutrino oscillations. 

(1) In the derivation of the oscillation probability in the framework of the wave packet 
approach, the integration of the squared modulus of the transition amplitude over time 
is involved. This is usually motivated by the fact that in most experiments the neutrino 
emission and detection times are not measured. One might therefore wonder if, when 
considering the experiments like K2K or MINOS in which the neutrino time of flight is 
actually measured, the integration over time is redundant or even illegitimate. 

In reality, the integration over time is necessary in any case due to the fact that the 
neutrino emission and detection processes are not instantaneous but rather take a finite 
time. This is reflected in the finite temporal extension of the neutrino wave packets at ~ 
a^^ ~ ax/vg (note that at depends on both the production and detection processes). One 
has to remember that even in the experiments like K2K or MINOS the neutrino production 

^^Note that in general this integral is not an L- independent constant, unlike the normalization integral 
(21). Still, it is Lorentz invariant. We encourage the reader to check that explicitly for Gaussian wave 
packets (see Appendix B). 



36 



and detection are not instantaneous but always take place within certain time windows. 
Consider, e.g. the K2K experiment [46]. The neutrinos were produced in the decays of 
pions originating from 1.1 /is proton spills, and the time window for the neutrino detection 
was 1.5 us. These time intervals are much shorter than the neutrino time of flight (~ 833 /is), 
but much longer than the time extension of the neutrino wave packets at ~ 10"^ fis}^ Thus, 
the time windows in the K2K experiments exceeded significantly the characteristic times of 
neutrino emission and detection, which justifies the infinite-limits integration over time. 

(2) It was pointed out in Sec. 3.1 and discussed in Sec. 3.2 that in general the mean 
momenta of the wave packets describing the neutrino production and detection states need 
not coincide. Pi Pi- This comes about because the conditions of neutrino production and 
detection are in general different. They are determined by the features of these processes 
(e.g., by the energy dependence of the production and detection amplitudes, the energy 
thresholds, the selection of certain energy intervals of the accompanying particles, etc.). 

At the same time, in the QFT-based approach to neutrino oscillations the neutrino is in 
the intermediate state and its evolution is described by a propagator characterized, in the 
momentum space, by a single momentum. The correspondence between the wave packet 
and QFT pictures can be established by noticing that in the wave packet formalism the 
functions Gi play to some extent the same role as the propagators in the QFT approach 
(see below), and that there are two different representations for these functions. These 
representations, given in eqs. (20) and (26), correspond to two different orders of integrations 
over the momenta of the production and detection states on the one hand and projection 
of the evolved produced state onto the detected state, which involves an integration over 
the spatial coordinate, on the other hand. In one case we first integrate over the momenta 
of the produced and detected neutrino states, thus forming the wave packets with certain 
mean momenta pj and p-, and then project the evolved produced state onto the detected 
state, as in eq. (20). The difference of the mean momenta then appears through the factor 
g^Pi-pDC^-i') in the integrand. This factor suppresses the projection if the average momenta 
are significantly different. In the other picture, which is illustrated by eq. (26), one considers 
the projection of the plane waves composing the production and detection wave packets and 
then integrates over the common momentum p. This common for the produced and detected 
states momentum (the same for and fP) appears as a result of the integration over the 
spatial coordinate x involved in the projection procedure, which leads to a delta function 
that eliminates one of the momentum integrations. The difference of the mean momenta 
Si = Pi — Pi appears then as a relative shift of the arguments of the shape factors in the 
momentum space. 

This second picture with a single integration over the momentum corresponds more 
closely to the one in the QFT approach. As we have already mentioned, the function Gi 

^^Thc length of the neutrino wave packet in the rest frame of the parent pion is essentially given by the 
pion lifetime: cr^ ~ ct-^ ~ 780 cm. In the laboratory frame it is about 1.5 orders of magnitude smaller due 
to the Lorentz contraction, which yields the above estimate for at ^ <Jx/vg (assuming that the time scales 
of the neutrino production and detection processes are of the same order of magnitude). 
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is a wave-packet analogue of the neutrino propagator of the QFT approach, which can be 
clearly seen from the representation of Gi in eq. (26). Unlike the usual propagator, however, 
it contains, in addition, certain information about the neutrino production and detection 
processes. 

(3) In this paper we have concentrated on neutrino oscillations in vacuum. In this 
case the wave packet approach (which implies a subdivision of the whole process into three 
stages: production, propagation and detection) and the QFT-based approach, which treats 
all three stages jointly, have comparable complexity, and the correspondence between them 
can be relatively easily established. In the case of neutrino oscillation in medium with a 
non-trivial density profile, the partition of the overall processes substantially simplifies the 
picture. As an example, solar neutrinos undergo flavor conversion inside the Sun, than 
propagate between the Sun and the Earth and then oscillate inside the Earth; the evolution 
of neutrinos has a rather complicated character, even when neutrinos are treated as being 
strictly real (on-shell) particles. In this case using the wave packet approach is essential, 
and it is important to understand its domain of applicability and limitations. 

(4) The coherent propagation condition (32) is the requirement that the wave packet 
separation upon propagation over the distance L be small compared to the effective spa- 
tial length of the wave packet ar^. It can be reformulated in the momentum space as the 
requirement that the variation of the oscillation phase within the wave packet due to the mo- 
mentum (or energy) spread be small. Indeed, since Acpst = (Am^/2p)L, the latter condition 
gives 

dA(l)st dp 
dp dE 

which is equivalent to (32). 

At the same time, the interaction coherence condition (33), which requires the variation 
of the oscillation phase along the wave packet in the coordinate space to be small, translates 
in the momentum space into the condition of no (or little) separation in the momentum 
variable between the wave packets describing different mass eigenstates. 

Mathematically, the loss of coherence due to the spatial separation of the wave packets 
is reflected in the lack of overlap of the factors Gi and Gk in the integrand of the coordinate- 
space integral representation (39) for Iik{L), and in fast oscillations of the integrand in the 
momentum-space integral representation (31). For the interaction coherence condition the 
situation is converse: its violation is reflected in fast oscillations of the integrand in the 
coordinate-space integral (39) and in the lack of overlap of the f^'^ and /^'^ factors in 
the integrand of the momentum-space integral (31). Thus, there is certain duahty between 
the wave packet descriptions of neutrino oscillations in the coordinate space and in the 
momentum space. 



cte = -TT^L — < 1 

Vg 



(95) 



38 



10 Discussion and summary 



In the present paper we discussed a number of subtle issues of the theory of neutrino 
oscillations which are still currently under debate or have not been sufficiently studied 
yet. For each issue, we were trying to present our analysis from different perspectives and 
obtained consistent results. We have developed a general wave packet approach to neutrino 
oscillations, which does not rely on the specific form of the wave packets. Our approach is 
formulated in terms of the generalized wave packets whose effective characteristics depend 
on the properties of the neutrino production and detection processes. In particular, we 
took into account that the momentum distribution functions of the produced and detected 
neutrino states are in general different, and explored physical consequences of this difference. 

We have also developed a new approach to calculating the oscillation probability in the 
wave packet picture, in which we changed the usual order of integration over the momenta 
(or energies) and summation over the mass eigenstate components of the wave packets 
representing the flavour neutrino states. This allowed a new insight into the question why 
the generally unjustifled "same energy" and "same momentum" assumptions lead to the 
correct result for the oscillation probability. 

We have also presented an alternative derivation of the equivalence between the results 
of the sharp-energy plane wave formalism and of the wave packet approach in the case of sta- 
tionary neutrino sources, as well as discussed the applicability conditions for the stationary 
source approximation. 

Below we give a short summary of our answers to the flrst seven questions listed in the 
Introduction. 

(1) The standard formula for the probability of non- averaged neutrino oscillations is ob- 
tained if the neutrino production and detection processes are coherent and the decoherence 
effects due to the wave packet separation are negligible. Under these conditions the addi- 
tional phase A0', which characterizes the phase variation along the neutrino wave packet, 
is negligible. The "same energy" and "same momentum" assumptions, which allow one 
to nullify this additional phase, are then unnecessary. These assumptions still lead to the 
correct result because their main effect is essentially just to remove this extra phase. 

An alternative explanation of the fact that the "same energy" assumption gives the 
correct result comes from the observation that in going from the oscillation probability to 
the observables such as event numbers, one has to integrate over the neutrino spectra. As 
discussed in [35, 14] and in Sec. 6, for stationary neutrino sources this is equivalent to integra- 
tion over the energy distribution within wave packets (provided that this energy distribution 
coincides with the spectrum of plane- wave neutrinos). In Sec. 4 we have shown that the 
integration over the spectrum of plane-wave neutrinos is just a calculational convention in 
the wave packet approach, which does not involve any additional approximations. 

(2) Quantum-mechanical uncertainty relations are at the heart of the phenomenon of 
neutrino oscillations. For neutrino production and detection to be coherent, the energy 
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and momentum uncertainties inherent in these processes must be large enough to prevent 
a determination of the neutrino's mass. These uncertainties are governed by the quantum- 
mechanical uncertainty relations, which also determine the size of the neutrino wave packets 
and therefore are pivotal for the issue of the coherence loss due to the wave packet separation. 

(3) The spatial size of the neutrino wave packets is always determined by their energy 
uncertainty aE- Note that aE is the effective uncertainty which depends on the energy 
uncertainties both at neutrino production and detection. 

(4) The coherence propagation condition ensures that the wave packets corresponding 
to the different mass eigenstates do not separate to such an extent that their effects can no 
longer interfere in the detector. This condition is therefore related to the spatial size of the 
neutrino wave packets, which is determined by the neutrino energy uncertainty aE- 

At the same time, the localization condition, which has to be satisfied in order for the 
neutrino production and detection processes to be coherent, is determined by the effective 
momentum uncertainty ap, which depends on the momentum uncertainties at production 
and detection. In the wave packet approach, the neutrino energy and momentum uncertain- 
ties are related to each other due to the on-shellness of the propagating neutrino, whereas 
in a more general quantum field theoretic framework they are in general unrelated. 

(5) The wave packet approach (or a superior quantum field theoretic one) are necessary 
for a consistent derivation of the expression for the oscillation probability. Once this has 
been done, wave packets can be forgotten in all situations except when the decoherence 
effects due to the wave packet separation or due to the lack of localization of the neutrino 
source or detector become important. Even in those cases, though, the decoherence effects 
can in most situations be reliably estimated basing on the standard oscillation formula and 
simple physical considerations. The wave packet approach is, however, useful for analyzing 
a number of subtle issues of the theory of neutrino oscillations. 

(6) The oscillation probability is independent of the production and detection processes 
provided the following conditions are satisfied: (i) neutrino emission and absorption are 
coherent, and decoherence effects due to the wave packet separation are negligible, and 
(ii) the energy release in the production and detection reactions is large compared to the 
neutrino mass (or compared to the mass differences). Note that if the condition opposite to 
(i) is realized, the probabilities of flavour transitions also take a universal form, as in that 
case they simply correspond to averaged oscillations. 

(7) The stationary source approximation is valid when the time- dependent features 
of the neutrino emission and absorption processes are either absent or irrelevant, so that 
one essentially deals with steady neutrino fluxes. Integration over the neutrino detection 
time then results in the equivalence of the oscillation picture to that in the "same energy" 
approximation. 

We have also demonstrated that, although calculations in the wave packet approach 
to neutrino oscillations are not manifestly covariant, they nevertheless result in a Lorentz 
invariant oscillation probability. This holds both in the case when the coherent produc- 
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tion/detection and coherent propagation conditions are satisfied, so that the oscillation 
probability takes its standard form, and in the general case, when these conditions are not 
necessarily fulfilled. 



• S D 

Appendix A: Dependence of f- {jp) and G^{L — Vg^t) on 
the index i 

For ultra-relativistic neutrinos to an extremely good approximation 

f-'''ip) = f'''ip,Em, (Al) 

i.e. the momentum distribution functions fi'^{p) depend on the index i only through their 
dependence on the neutrino energy Ei{p) = >y p'^ + mf . The same applies to the shape fac- 
tors Gi{L — Vgit). Indeed, these functions depend on i through the neutrino mass dependence 
of the phase space volumes and of the amplitudes of the neutrino emission and detection 
processes. The phase space volumes depend on rrii only through Ei{p); the production and 
detection amplitudes can in principle depend directly on rrii, e.g., due to the chiral suppres- 
sion, as in the case of tt^ or decays. However, the corresponding contributions to the 
total amplitudes are completely negligible compared to the main contributions, which in 
this case are proportional to the masses of charged leptons. 

Since Vgi = Ei{py — m'f/Ei{p), eq. (Al), in particular, means that in the limit when 
r = 1 (i.e. Vgi = Vgk) one has fi'^{p) = fk^{p) and Gi = Gk- 

Appendix B: Lorentz invariance of the oscillation prob- 
ability and Gaussian wave packets 

We illustrate here how our general proof of Lorentz invariance of the oscillation probability, 
given in Sec. 8, works in the particular case of Gaussian wave packets. To this end, we shall 
prove the invariance of the integral 

hk = J dt Gi{L - Vg,t)Gl{L - vgkt) . (Bl) 

Lorentz invariance of the oscillation probability (22) will then follow from the definition of 
lik in eq. (23) and Lorentz invariance of the phase difference A0 (24). 

The normalized Gaussian wave packets of the neutrino mass eigenstates in the laboratory 
frame are [25] 

^1/2 
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Here we have taken into account that the lengths a^i of the wave packets corresponding to 
different mass eigenstates z/j are in general different. Even though the dependence of a^i 
on i is expected to be extremely weak, we take it into account here, to be in line with the 
transformation law (86) which depends on the group velocity of the wave packet (one can 
neglect this dependence at the end of the calculation). Note that the exponential factor in 
(B2) is Lorentz invariant (cf. eqs. (94) and (86)), whereas the pre-exponential factor is not. 

Consider the integral lik with Gaussian wave packets Gi{L — Vgit) from (B2): 
Direct calculation of the integral yields 



L'2 (Vgi-Vgkf 



Hk 



(^l^gk + (^Ik^li 



-'gi'^gk \ 4(0-^ .d'-^, +0-2, v^.) fi„^„'^ 

■J ^ 1 r) ^ XI gk xk gi^ p OCTrrUn 



3 



(B4) 



Here in the last equality we took into account that for relativistic or quasi-degenerate neu- 
trinos one can set = (Jxk = and also neglect the difference between Vgi and Vgk in the 
pre-exponential factor, while retaining the difference Avg = Vgi — Vgk to leading order in the 
exponent, where it is multiplied by the macroscopic distance L. The resulting expression 
can be written as I^k — e"^'^/'™'^-' , i.e. lik is the overlap integral, which takes into account 
the suppression of neutrino oscillations when the baseline exceeds the coherence length and 
the wave packets corresponding to different mass eigenstates separate and cease to overlap. 

Let us now discuss the properties of lik with respect to the Lorentz boosts. By changing 
the integration variable, eq. (B3) can be rewritten as 

^ [V9^(7xk) / dS_ ^ 3^ ^ 

{27rvgkaxiy/^ J a^k 

where s = L — Vgkt. As follows from eqs. (94) and (86), the integration measure ds/a^k is 
Lorentz invariant, and so is the exponential factor in the integrand (note that it coincides 
with that in eq. (B3) and has just been rewritten in terms of the new integration variable); 
therefore, the integral in (B5) is Lorentz invariant. In a moving frame, the normalization 
factor of the Gaussian wave packet (B2) has to be re-calculated (see the discussion below), 
which gives axi^k — 0"xj,fc5 ^si.fc ~^ ('^gi,k ~ faking this into account and making use of 
eqs. (85) and (86), it is easy to show that the factor in front of the integral in eq. (B5) is 
also invariant under Lorentz transformations. Thus, the quantity I^k is Lorentz invariant. 

It is instructive to check explicitly how Lorentz invariance works in the case of Gaussian 
wave packets. Consider the overlap integral /^'^ in a reference frame where the neutrino 
source and detector are boosted with a speed u in the direction of motion of the neutrino. 
It should be stressed that one cannot obtain I-j^ by simply replacing the unprimed quantities 
in eq. (B4) by the primed ones - the integral has to be calculated anew. This is related to the 
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important difference between the properties of L and V . These are the distances neutrinos 
propagate, respectively, in the laboratory frame and in the moving frame, and they are 
defined as the differences between the corresponding coordinates of the neutrino absorption 
and emission points. However, while L is a fixed distance between the neutrino source and 
detector in the laboratory frame, which does not depend on the neutrino propagation time 
t, the corresponding distance L' in the moving frame is not fixed and actually depends on t\ 
as follows from eq. (88). In particular, L' is not just the usual Lorentz-contracted value of 
L (equal to Lj^^; the latter would correspond to the difference of the detector and source 
coordinates in the moving frame at the same time t', whereas L', by the construction of the 
effective wave packets, is the difference of the detector coordinate at the time t' and the 
source coordinate at t' = 0. 

The dependence L'{t') can be found by substituting t from the second equality in (88) 
into the first one, which gives L'{t') = L/'-fu + ut' . Taking this into account and replacing 
in eq. (B3) all unprimed quantities by the primed ones gives, after a simple calculation, 

(Note that the normalization factor for G[{L' — v'^^t') also had to be re-calculated, since the 
normalization of Gi{L — Vgit) in (B2) was obtained for t-independent L). From eqs. (85) 
and (86) one finds 



'^l{l+UVgi){l+UVgk) 



I I '^gi ^gk 



'Vgi ^gk 



7^(1 + UVgi){l + UVgk) 



12/ I \2 , l2 / I \2 ^ xi'^ gk ~ ^ xk'^ gi /-r,^\ 

Substituting this into (B6) and comparing the result with (B4), one can readily make sure 
that /-fc = iik. 
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